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Abstract. We prove the exact asymptotic 1 — + o(l)) j\pn 

for the probability that the underlying graph of a random mapping of 
n elements possesses a unique highest tree. The property of having a 
unique highest tree plays a crucial role in the solution of the famous 
Road Colouring Problem [S] as well as its generalization in the proof of 
the author’s result about the probability of being synchronizable for a 
random automaton [T]. 


1 Preliminaries 

Given an integer n > 0, by the probability space of random mappings of n 
elements we mean the set of all n" mappings from Q = {1, 2,..., n} to itself with 
the uniform probability distribution. Each mapping / S Sn can be represented 
via the directed graph g{f) with constant outdegree 1. The graph g{f) has the 
vertex set Q and the edge set E = {(p, f{p)) \ p S Q}. Since this is a one-to-one 
correspondence, we identify with the probability space of random digraphs 
with n vertices and constant outdegree 1. 

Each digraph having common outdegree 1 consists of cycles and trees rooted 
at these cycles. Let T be a highest tree of g and h be the height of a second by 
height tree. Let us call the crown of g the (probably empty) forest consisting of 
all vertices of height at least h+1 in T. A digraph g on Figure [1] has the unique 
highest tree rooted at state 2. The height of a second by height tree is 2 whence 
the crown of g is the forest consisting of the states 7,8,9,14 of height al least 3. 
Denote p = ^ -~ X ~ results are as follows. 

Theorem 1. Let g G Sn be a random digraph and H be the crown of g having 
r roots. Then |i7| > 2r > 0 with probability 1 — 0{l/y^), in particular, a 
highest tree is unique and higher than all other trees of g by 2 with probability 
1 - 0{l/^/n). 

Theorem 2. The probability for random mapping on n elements of having ex¬ 
actly two highest trees is -^(1 + o(l)). 

Corollary 1. The probability that the underlying graph of a random mapping 
of n states has a unique highest tree is 1 — -^(1 + o(l)). 


In order to apply aforementioned results in [I] , we need an easy generalization 
of Theorem[T] Given a digraph g G En and an integer c > 0, let us call a c-branch 
of g any sub-tree of a tree of g with the root of height c in g. For instance, the 
original trees are 0-branches. Let T be a highest c-branch of g and h be the 
height of the second by height c-branch. Let us call the c-crown of g the forest 
consisting of all the vertices of height at least h + 1 in T. For example, the 
digraph g presented on Figure [1] has two highest 1-branches rooted in states 
6,12. Without the state 14, the digraph g would have the unique highest 1- 
branch having the state 8 as its 1-crown. 

Theorem 3. Let g G be a random digraph, c > 0, and H be the c-crown of 
g having r roots. Then \H\ > 2r > 0 with probability 1 — in particular, 

a highest c-branch is unique and higher than all other c-branches of g by 2 with 
probability 1 — 0{l/^/n). 



Fig. 1. A digraph with a one cycle and a unique highest tree. 


We start from the proof of Theorem [5] and obtain the rest as consequences 
of the proof. To prove Theorem [2] we use the following scheme. The probability 
distributions (p.d.) on random digraphs from En having N cycle vertices can 
be considered as the p.d. on the set of random forests Fn,N with N roots and 
n vertices averaged by the p.d. of the number N of the cyclic states of random 
mappings. In its turn p.d. on random forests can be considered as p.d. on critical 
Galton-Watson branching processes under a condition on the total number of 
particles in the process. This idea of using Galton-Watson branching processes 
plays a crucial role in probability analysis of such combinatorial objects like 
random trees, forests and mappings. For an introduction to this theory, we refer 
the reader to [H Section 2]. 

Unfortunately, in such a reduction there is an obstacle of using the theory 
of Galton-Watson branching processes out of the main range, namely, formulas 
become too rough when the number of cycle vertices is “big” and the height of 
random forests is “small” with respect to n. 

Our proof is organized as follows. First, in Section [2] we represent p.d. on 
mappings via p.d. on forests averaged by the number of cycle vertices N, and then 
reduce the range for N. Next, in Section [3] we use the standard representation 
of p.d. on random forests via branching processes. Using this representation and 
an independent approach, in Section |4] we reduce the range for the height of the 


forests with respect to both n and N. Then in Section [5] we present a theorem 
concerning branching processes for our case in the main range of parameters, 
and based on this theorem and results from the previous sections, we prove 
Theorem O Finally, in Section [5] we get remained results as consequences of the 
Theorem [2] proof. 


2 Number of the Cyclic Vertices 


Given a digraph g € En, denote by X{g) the number of the cycle vertices and by 
ft{g) € Fn-\(g)^\(g) its forest, that is, the union of the trees of g. Given d > 1, 
denote by the set of digraphs g € En having exactly d+1 highest trees, by 
Bn the set of digraphs having at least two highest trees and by the set of 
digraphs which do not satisfy Theorem [T] Given 0 < t < n, d > 0, denote by 
B'^ ( the subset of mappings from B'^ whose second by height tree has height t. 
Given 0 < N < n, denote by t — Fn,N the set of forests of digraphs from 
Bn t with N cycle vertices. Due to the definitions, we have 


n — 1 

Since the definition od depends only on the forest of a digraph, we also have 

n 

(Bit) = E (Hg) = N)P{Bt^^^,). (1) 

N=0 


The following formula is well known (see e.g., [3l Lemma 3, Section 3]). 


P{X{n) = N) 


N{n-1)\ 
n^{n — N)\ 


( 2 ) 


The following corollary along with o allows us to reduce the range we have to 
consider for the cyclic vertices number. 

Lemma 1 (see Sec. [3 for the proof). 


P{-^ < A(n) < 4\/ n In n) = 1 — o{l/y/n). 


The following corollary easily follows from the proof of Lemma [T] 

Corollary 2. 

i^.ei:„(A(n)=iV) = ^^(1 + 0(1)), 
y/n 

where z{N,n) = N/y/n varies in the range 4\/lnn). 





3 Random Forests via Branching Processes 


We use the standard representation of random forests distribution via critical 
Galton-Watson branching processes following [3] . Denote by the probability 
space of all critical Galton- Watson branching processes with the offspring p.d. 
(e“^/r!)j.>o starting with N founding ancestors (or, equivalently, consisting of 
N independent branching processes p} ,iJ?,... from G^) where each particle 
independently has probability e~^/r\ of producing r offsprings in the next gen¬ 
eration. The idea of such representation is that after a branching process with N 
founding ancestors stabilizes, we obtain a forest having N trees corresponding 
to the processes and the induced p.d. on the set of random forests 

is uniform (see [3l Section 2] for details). 

Following [5], denote by F{z) the generating function for the offspring p.d. 


+ 00 


^W = E;h = 


.Z-1 


r—0 


( 3 ) 


Given /i G G^ and t > 0, fj,{t) (resp. i^(t)) denote the number of particles 
in the process /i in generation t (resp. less than t) and by v denote the number 
of particles in p, after the extinction moment of /r, that is, v = vir^p)), where 
r(/i) is the index of the first empty generation of p. Notice that the height of 
the induced forest exceeds the extinction moment of /r by one. 

For processes with multiple founding ancestors, the subscript is used to de¬ 
note the number of founding ancestors (0-generation particles) and the super¬ 
script is used to denote the index of the ancestor in 0 generation. In particular, 
UN G G^ denotes the total number of particles in Galton-Watson process start¬ 
ing from N founding ancestors and iz* denotes the number of particles of the 
i-th founding ancestor. By the definitions. 


N N 

UN = '^ u\ UN,t = UN{t) = ^ u\t) 

i=l i=l 


where all F are independent. 

Now we determine the conditions on branchi:M processes corresponding to 
the sets of digraphs B!^ To simplify notations 111, suppose that .. are 

ordered by descending of extinction moments r^, r^,..., r. Given iV > 0, t > 0, 
denote by ^ the set of processes from G^ such that 

fj.^{t) >0, — u^(t) < fJ-^{t), t -I- 1, fi^{t -I- 1) = 0 for j > 3 (4) 

and for d > 0 by Af^ ^ ^ the set of processes from G^ such that 
H^{t) > 0, fi^{t-\-r) =0, T* = t-|-l, t for 2 < * < d-|-l,j > d-|-l. (5) 

Also denote Af^ ^ = and AN,t = ^d=i^N,f 


^ For calculations, we cannot assume this and the order matters. 



For d G {_, 0,1,..., — 1}, due to the definitions of ^ and Aj^ the 

standard reduction to branching processes yields 


P{Bn,N,t) — P{^N,t \ i^N — n + N) — 


— n + N) 


P{vn =n + N) 

whenever P{vn = n + N) > 0. The following lemma is given in [3]. 


( 6 ) 


Lemma 2 (Lemma 6, Section 3 [3]). P{yN = n + N) = (1 + o(l)), 

for 0 < zo < N/^/n < Zi. 


The proof is based on the following equation 

N 


P{vn =n + N) = 


{n + N)^j2'iTB{n + N) 


e-2B(„+N)(l + o(l)), (7) 


where i? = _F"(1) = 1 in our case (see ([3])). In its turn, this equation is based on 
O Theorem 2, Section 1] which holds for any n,N ^ +oo. Hence from ([7l) for 
N = o{n), we have the following corollary. 


Corollary 3. P^un = n + N) = ^=^(1 + o(l)), 

where z{N,n) = N/y/n varies in the range (n“°'^®, 4-\/lnn). 

By Corollary [31 we have 

, P(h'N = n:A%f) j 

Pf^F^AKN,t) = ^ = (1 + o{l))n^P{v^ = n; (8) 

for z = N/y/n G 4\/ln n) and by Lemma [H we have 


4\/ n In n 

Pgei:jB^,t)= E P(Mn) = N)P(B%^,)e^(n), 

N= ^ 

where ei(n) = (1 + o{l/y/n)) 

Putting together ([9]),(|8]) and Corollary |2l we have 


(9) 


4 a/ n In n / n, 

n /T^d y V- £i{n)ze 2 

Pgf^sABn^t) - 2^ - 7= - 

___ V ^ 


4V n In n 

= ei{n)y/n E = n; A% t). 

N= ^ 


Thus in order to prove Theorem |2l it remains to show that 


n 4\/n In n 

E E = n■,A]yf^) = ^{1 + o{l)). 

t=l N=^ 


( 10 ) 















4 Height of Random Forests 

Now we want to reduce the range we have to consider for the parameter t. For 
this purpose, we use the following basic equation (see e.g m 

P(/ri(t)=0) = l-^(1 + 0(1)). (11) 


Lemma 3. If t > 0,N > 2tlnn, then P{vn = = o(;^)- 

Proof. Due to the definition of we have 

P{i>N = n-,Ais!,t) < P{AN.t) < = 0 )^“^. 

It follows from (HU that there exists to > 0 such that P{fj,i{t) = 0) < 1 — 
for t > to- Denote po = minjxtg P{pi{j) = 0). Clearly po < 1- Hence for t < to, 
the lemma trivially holds. For t > to, we have 

1 QQ 1 

P{,^M = n-, AM,t) < 0(1)(1 - —< 0(l)e-3-8i"" < o( —). 

t U'^ 

Following [7] denote by Tn^u the number of mappings of n element set of 
height at most h, by Lj(x) the j-th iteration of the function xe^ and by pj the 
unique real positive solution of the equation Lj(x) = 1. 

Corollary 4 ([Tj). El + pi + pip 2 ^ -h P 1 P 2 ■ ■ 

In another work, Grusho [H Lemma 4] found the following asymptotic for 
where m/n^ —>■ C > 0 for 0 < /3 < 



Lemma 4. For each e > 0, PgeSr^iBn,t) = O(^). 

Proof. It follows from Lemma [T] and Lemma [3] that 

^/n/{21nn) n .. 

E E = o(^). 

N^l ^ 

Using that pm is monotonic descending, for < t < due to m and 

Corollary ID for m = and n big enough (such that t + r< m), we have 


P{Bn,t) < 


Tn,t+r ^ {^ + f + 0{^))) ^ 

- y ) gn - 



2 


The Corollary 0] is given for constant h but it is not used anyhow in the proof. 







5 The Main Range 


We write (^i(n) ^ (1)2{n) if = (1 + o{l))(j) 2 {n)- The following theorem is 

based on an analysis of presented in Section [71 

Theorem 4. Denote P = P{9,n) = \J—^iOjn; Then for N G (-^,4\/‘nlnn),t G 
n),iV/t < 2Inn, we have 


P(vm — n\ A^f f) 


AkU I n In* n (1 — 

■ {^ + ^ 3 {t:—))d9 + o{ — ), (13) 
n 


where |e 3 (t,M)| < 


0(1) (lnn(v1^+ ^)) 


2 


and for any constant r > 0 


P{vN = n; Ajv,t,r) < 0(l)P{vN = n; (14) 


Proof. Let f{z),ft{z) be the probability generating functions (pgf) of iz and 
V, /i(t) = 0 resp., that is, 


+ 00 


+00 


Mz) = = 0)z", f{z) = Y.P^^ = 


n—0 


n—0 


First let us make the following remarks. 

Remark 1. For given original indices of the pgf of [vn] Aiq^t,r) (and 

P[vn = n;Alf^)) is equal to fl^~'^{z)(l)i{z)(l) 2 iz), where 4 ‘ 1 t ^2 are the pgfs of 
{v'^\AN,t,r),{v'^]AN,t.r) (for r = 1) resp. 

Remark 2. For any r > 1 there are at most N{N — 1) choices for the indices of 
and for fj, G Aj^ t.r and there are exactly choices for and for 

T e 

Notice that if g{z) is the generating function of a random variable X taking 
positive integer values, then y(0) = ^(e*®) is the characteristic function of X. 
Due to the inversion formula, we have 


Pirn = n; ^ £ e-^®"/f-"(e^')(/t+i(e^®) - /t(e*^))^d0 ^ 

]\t 2 pTTn 

- ^ (e^" )(/t+i(e*" ) - ft{P-)fd9. (15) 

47777, J _7rn 

Notice that 

|/t+r(e*" ) - /t(e*")| < P{t < T{fi) <t + r)< 


( 16 ) 











We consider three ranges for 6 . Due to LemmafTUl there are positive constants 
e, ^ 2 , Cl, C 2 such that 


|/t(e*")| < 1 - i 52 \/|-| + cit 
'' n 


— 1—C2t 


(17) 


for \0\ < e. Denote A = 2^. Then > 2 cit-ie-^"V^ for |0| > A and 

n big enough. Hence using (HZl and m, we have that 




< 

N 


< 




/ en / I f) \ 1 N /~ 

h_ 0.5^2^1-11 de<^J^ (18) 


where C 3 > 0 is some constant. Notice that 


I 2 

e-<iVSd 0 = [y = q^^]= / ve-^dv = 

r Jg^c 


2 {u + l)e .q^ 

Qr2 Ij-v/d’ 


Substituting q = c = A, d = en, we get that (IT^ is upper bounded by 


0 ( 1 ) / \/n'/A 


Va ^ 0 (l)n 


— e ^ < 

y N m) - t4 


C 3 In n _ 


= «(Z4)- (19) 


Due to Lemma 1 there are constants gi < 1,(72 < 1 such that |/t(e*“)| < 
qi+ q\ for e < |u| < tt. Hence we have that 



(e*-)(/t+r(e^") - ft{e^-)?de 


<^{di + <lir <0{\)^,q^ ( 20 ) 


for some constant (73 S { 1 , 02 ). 

Notice that —>■ 0, for \9\ < A,t G (n°'49,n). The (fT51) follows from 

(fT^ . (fT^ . (l20l) and Corollary [7| (IT4l) follows from the proof and Remarks [U [2] 


Due to Theorem |4j it remains to compute the sum for all integers t, N such 
that tG (n°-^®,n),iV G TOin(4\/nlnn, 2t Inn) of 
















First let us sum it up by N. Since ^(iV) is smooth and positive for > 0, 
we have that 


;tlnrt |^2tlnn 2tlnn / pN+l \ 

= / (t^{z)dz + ^ (i){N) - / (j){z)dz = 

N=^ J 

/ 2tlnra 2tlnn 

(l){z)dz+ ^ (j)'{N + S'j^) = 

N=^ 

p2tlnn p2tlnn 2t\nn 

= / (j){z)dz + / (j)'{z)dz + ^ (j)"{Ns) = 


N= ^ 


f j.2t\nn \ 2t\nn 

/ (j){z)dz + (j){2t\nn) — (j)"{Ns)^ 

J N=^ 

„2tlnra 2tlnn 

- / (l>{z)dz+ (t>''{Ns), (21) 

J ^ 


N— ^/n 


where S'jy G (0,1),TV < iV^ < + 1. 

Denote u(0) = t^J^,a{6) = jd following lemma follows from 

[3] [Lemma 5, Section 3] but we reproduce it here for completeness. 

Lemma 5. For some positive constant c, Re{a(9)) > 

Proof. One can check that 

^Re{a) _ l-e-2“ + 2e-" 

u 1 — 2e““cos(u) + 

It follows from here that (/)(m) = as the function of u is positive, continuous 

and lim„_>o 4>iu) = +oo and limu_,.o (f>iu) = 1. The lemma follows. 


From the other hand, we have that 


^(^)l < 


i(l 


t(l - e-“) 

The latter inequality follows from 


< 


3(m + 1) 


1 < 


1 — e 

for each a; > 0. Since for |0| < ” , 


|e3(A-))l < 0(1) Inn 

n \ 


— < 3(a; + 1) 


( 22 ) 


( 23 ) 


6 , 1 
~\^~y ~n 

n vT, 


<o(i), 


2 











for the second term of m, we have 


2t In n 

I E 

2t\nn 

i: 

N=^ ~~P~ 

2tlnrt - " " 

2t In n n In^ n 

i: /:... 

N=^ ~~t^ 


< 0(1) In^ n E 

N=^'‘ 


'4p-2t/3 




(_(27V, _ N!a)a + 2 - 2fV,a) e-0)4 


-Nsa 


dO < 


(Nsa - 2f 


o4 -2t/3 

P ^ ^-N,a 

^4g-2t/3 


(lnn(. + l) + 2)^ ^EVo)4 ^~^^“ 


d9 < [^,N < 2t\nn] 
dO < [Lemma [5] 


n In^ n 

(u + l)^M‘*e 


Jo 

t4(l - e-“)4 

- ^(l-e-“) 


2tlnn . " In** " , , ,x6 

<0(l)ln2n E [ ‘ ^^^e-^“e-"^d0<[e-"^<l] 


Ar= -J/n " 


< 0(l)ln3n 


(u+1)® _ 2 u,a tjo 2nw 
-- e d0 = [dS = —^du] = 


Ard 

= 0{l)\n^n / 

Jo 








-■ (24) 


Summing up 0(4)”'^ " by t S (n°-^®,n), we get that it is bounded by 0 (^;t^). 
For k > 0 and a having the positive real part, we have 


t. 1 t efe 

x^e-^^dx- ' k 


-y/c+1 


’^dx = 


fc+i ■ 


(25) 


where is a constant which depends on k and for k = 2 


^+oo 


c^e-^^dx = 


_™ x‘^ + 2x 2 . 


^fc+l 1+00 ^fe + 1 ■ 


(26) 


































For the first term of (1^ . by (051) . we have 


2ilnn /*+oo 

(j){z)dz ^ / (j){z)dz 
0 

71 In ' 


^ /34e-2‘/5 


I A o n -t/ 3^4 3 ^^ ^^(t,—))d 0 — 

I n In* n ZTTfl — 6 Tl 


I In^ 7 


-i9f 




_*e/3e 2*/3(i + e ^ 

u„4„" 2^(l-e-‘/3) 


,0/3e 2‘/3(i + e */3)3 0 

2^(l-e-*/3) "3(i, (27) 


Recall that m( 0) = ty-^- We have 


|e3(t,^)| < 0(1) hnn < 0{l)^-^{\u\ + 1)^ 


Using that Re{t(3) = u, we can upper bound the second term of (l?7l) as follows. 


0 ( 1 ) 


/3e-20(i + e-0)3 ff 

(1-e-O) ^3(i,^) 


d0 < 


<o(i) 


ue ^“(1 + e “)^ In^ n 


t(l - e-“) t 


(u + 1)^ 


d0 < [dd = 


An^ 

< 0{n In^ n) / 

Jo 


y2g-2«(i _^g-«)3 


f^(l — e“") 


■(« + l)^ 


du < [— -r < U + 1] < 

“ yi - e-“) “ ^ “ 


< + l)"e-^“(l + e-'^fdu < 0(^^), (28) 


The sum of "• ^ by t G (n°-^®,n) is upper bounded by 


































Let us now consider the first term of (l?7l) for positive 9. 




27r(l — e 

(1 - z)ue-2“(i-d(l + 


n 


^In n 


2'Kt{l — e *)) 
i (1 - i)u2e-2“(i-d(l + f,-n(i-i)Y 
7r(l — 


dO = [d9 = —^du] 


du = Tpit). (29) 


Using the same trick as in (|2T]) for have 


n .n n 

^ / 'lp{t)dt+ ^ 

„n.4Q ^ „n.4Q 






(30) 


where ts G {t,t + 1). 

For the second term of (ISO)) , we have 






<0{l)n 


n / 

Yi 

nu\ 

. nu 

1 

1 1 



^ '¥ 



(1 - i)u2e-2“(i-d(l + e-“(i-0)3 


(1 - e-“(i-d) 


du < 


<o(i) 


n In^ n 




E 

t=n0-*3 • 


< 0(1) r 

Jn<^ 


< 0 ( 1 ) 


'•1" " u2 


^ + -^ + 1 ) m(m + 1) 


^du < 


n In^ n f n 

% + -r + 1 Mi < 


f5 t4 

n M n f n n 


+ TT + 1 


t4 Vi2 ' ^4 


< o(n-‘''y^). (31) 















Finally, for the first term of (1501) . we have 


/ ^it)dt = [v =- — ,dv =—] = 

J„0.49 r 


Ard n 


dv 


(1 - f)M2e-2“(i-*)(l + e-“(i-*))3 


27r(l - 

fin^n 

3 27rf(l — 

^ - e-^- )ge~^^{l + e“^)^ 




1 — i 


0 

^(1 —i) In^ n 


47r 


, In^ n 


1 — i 
Att 


+ 00 

E 

k=0 


l + 0{ -)-e 

n 


1 3 

+ 


27r(i + 1)(1 — e s) 

(i-»)an°' 


du = [g = (1 —i)^] 
dg = {A) 


+ 00 


ge-2s(l + e-9)3^e-'=®d5= (i?) 

1 


fe =0 


(fc + 2)2 (fc + 3)2 (fc + 4)2 (fc + 5)^ 


(1 + 0 ( 1 )) 


8 - 

7 

4 


1 

( 47r2 

827\ 

1 6 

4 

9 

16/ 

47r ' 

1 3 

144/ 


(32) 


where (A) follows from ^ = J2t=o d'" and = 1 + 0(12;;;^) for I 5 I < 2 In^ n 

and (B) follows from the equality ge~^^ = "). 

One can easily obtain for negative 9 the same expression with the factor 1 + i 
instead of 1 — i. Thus we get that 


EE P{vn — n] A]y j) — —(1 + 0(1)). 


i=l N=1 

Theorem [5] now follows from (1551) . 


(33) 


6 Conclusions 

Notice that Theorem [T] corresponds to digraphs i?„ which correspond to branch¬ 
ing processes Ajv,t- Hence, the lower bound of Theorem [T] follows from the fact 
that A]y t implies Aw,*- To prove the upper bound, we use inequality (1141) of 
Theorem |4] and the following lemma. 

Lemma 6 (see Sec. 0 for the proof). There exists some constant r > 0, 
such that for each t > ® the following inequality holds for n big 

enough. 

Pfu = rn;iy — u{t) < g{t)) < 

< 0 (l/n°' 2 i'i") + c 2 P{v = m(l + o{l))■^l{t) > 0 ; g{t + r) = 0 ). 




















The above lemma allows to neglect the difference in the definition of and 
J^N,t,r for proving the upper bound. The only thing we have to notice to use 
Lemma ini is that P{i'n = 'n]Aiq^t,r) is continuous by n in the main range of 
parameters, i.e. P(yN = ''^■(1 + o(l)); ^ P{vn = 'n;Apj^t,r) (this easily 

follows from the proof). The precision of Lemma |6] is enough due to the reduction 
([HI) and Corollary m 

In order to prove Corollary [1] consider the set of processes in having at 

/N\ f 3^G —\ ^ 

least 3 highest trees. Then, we would get a factor ( (i_e -*/^)2 ) instead of 

(h ( (i-e-tg) 2 ) ■ Since | \ < 0(l/t^), we would have the additional 

factor of order |^. Due to Lemma[31 we consider N < 2tlnn whence ^ = o(l). 
Corollary [T] now follows from Theorem [5] 

Since the generating function of /i 7 v(l) is given by 




+00 „■ 

iV —1^ = p^c-1) = 

pn\> ^ j\pN ’ 

3=0 3=0 

we get that P(|p,Ar(l)/iV — 1| > 1/2) = o{l/N^). Due to Lemma[I]this means 
that whp /rAr(l) is equal to N up to the constant factor, and the Theorem [T] 
holds for 1-branches instead of the trees. By induction, one can easily generalize 
these arguments for any constant c > 0 and Theorem [3] follows. 

7 Technical lemmas 


Lemma [TlP(-^ < A(n) < 4-\/n In n) = 1 — o[l/^/n). 

Proof. For N = n, P{X{n) = N) = ^ = 0(i). For n —>■ -|-oo,A^ < n, by 
Stirling’s formula fc! = (|)^-\/27rfc(l -I- o(l)), we have 

P(X(n) = N)= = (1 + o{l))N{{n - l)/e)"-^V^ ^ 

n^{n — N)\ {fri — iV)/e)"“'^-\/27r(n — N) 

NrP-^(l - i)"-i 


e^-in"(l - y n-N 

il + o{l))N 


n 


(l + o(l))x 


(34) 


e^n(l - ^)"“^ y n - N pn{x+(i-x)\’a(i-x))_ j.’ 

where x = ^. Denote (/)(x) = cc -|- (1 — a;) In (1 — cc) for 0 < x < 1. Using Taylor’s 
expansion In (1 — x) = —x — ^ — 0{x^), we get that 




(j){x) = x + {x - l)(x -I- + O(x^) = ^ + o(x^). 


(35) 


Since (/'(x) = 1 — 1 — In 1 — x > 0, (j){x) is growing. Hence by (l35l) for x > 4 



















Since also 


/1—X 


< y/n for 1 < < n, by (IMl) and (1551) . we get that 


P{X{n) > Ay/n Inn) < ny/ne = o( —). 


For N < Ay/n Inn, x = o{l) whence by (l34l) and ([35]), we have 
P(A(n) = TV) = 


(37) 


From here, we have also that 


P(A(n) <^)< ^ (1 + 0(1))-^ = o(^). 

n y/n 


The lemma follows. 

Lemma|6[ There exists some constant r > 0, such that for each t > m > 
® the following inequality holds for n big enough. 

Pfu = rn;iy — u{f) < g.{t)) < 

< 0(l/n°-^''"") + C 2 P{n = to( 1 + o(l)); ^(t) > 0; fj,{t + r) = 0). 

Proof. Since the distributions of vft) and /i(t + r) are determined by qi{t), for 
the left hand side, we have 

P{v{t) = m;v — v(t) < fi{t)) = 

+00 

= ^ = k)P{u{t) = m I fi{t) = k)P{u — u{t) < fi{t) I fi{t) = k) = 

k^l 

+ 00 

= ^ P{f^{t) = k)P{'ipt,k = m)P{uk < 2k), (38) 

k^l 

where ipt.k = v{t) \ qi{t) = k. For the right hand side, we have 
= m'; /r(t) > 0; ^{t + r) = 0) = 

+ 00 

= P{fj,{t) = k)P{i' = m I = k] = 0)P{fi{t-\-r) = 0 | fi{t) = k) = 

k^l 

+00 

= ^ P{t^{t) = k)P{ipt,k + Vk,T = m')P{qi{r) = 0)'". (39) 

fe=i 

Due to Corollary 121 P{vk < 2k) < C 2 e“^/® and due to (fTTl) . P{fJ,{r) = 0) > 
gi/5 for some constant r > 0. Thus it remains to show that nk,r = o{m) with high 
probability. Suppose that i'k,r > kD^ for some D > 0. Then at some generation 







i < r one of the particles in z-th generation must have at least D offsprings. 
Hence 

k 

P{vk,r > kD") < 

eD\ 

For k > In^ n, we have that 

P{vk < 2fc) < Cae"'""’"/® = 0(l/n°'2 

For k < In^ n, we have that 

P{i’t,k + Vk,r = m!) = P{i’t,k + Vk,r = m' I Vk,r < k'nP-^)P{vk,r < fcn°-®)+ 

+ P{'>pt,k + Vk,r = m' I iyk,r > kvP-^)P{vk,r > kn°-^) < 

/ , / / NNN In^ n 

< P{iPt.k = m(l + o(l))) + = 

= P{'>Pt,k = m{l + o(l))) + 

So, we have 

=m-v- v{t) < < 0(l/n°'2'''") + 

In^ n 

+ ^ P{li{t) = k)P{ilJt,k = rn)P{yk < 2k) < 

k=l 

< + o(n-°'^") + C 2 Piiy = m(l + o(l)); Az(t) > 0; + r) = 0). 

The lemma follows. 

Recall that f{z),ft{z) are the probability generating functions (pgf) of iz and 
= 0 resp., that is, 

+00 +00 

= 0)^", f{z) = '^p{v = n)z^ 

n—0 n—0 

and denote At{z) = f{z) — ftiz). The following lemma is proved in 
Lemma 7 (Lemma 1 [5]). If t ^ +oo,u —^ 0 such that ut —>■ 0, then 

+ (40) 

where a{u) = \/—2m (the branch with a positive real part is always chosen) and 
\e{t, u)| < clnt('ut + 1/t). 

Corollary 5 (Corollary 1, Section 2.3, Page 127 [3]L For each e > 0 there 
exists constant ga < 1 such that |Z\t(e™)| < ql for e < |u| < tt. 

Corollary 6 (Corollary 2, Section 2.3, Page 128 [^). There exists e > 
0, Cl > 0,ca > 0 such that |Z\t(e*")| < for |u| < e. 






It is proved in that for some c(u) such that |c(it)| < C 3 , 

f(e^^) = l — a(u) + c(u)a^(u), (41) 

Using Lemma [3 and Corollary [SJ by (HD) we get the following lemma. 
Lemma 8 . If t ^ + 00 , u —>■ 0 such that ut —>■ 0, then 


f^[en=[l-a{u) + c{u)a\u)--^-L^^{l + e{t,u))^ = 

( 1 _|_ \ 

= {\ +Ne:i{t,u))exp\^Na{u)-^-—^^^^y 

where |e 3 (t,u)| < 0(1) (|u| + \A{t,u)\e{t,u)) < 0(1) Int (|u| + . 

Lemma 9. For each e > 0 there exists constants gi < 1,(72 < 1 such that 
|/t(e™)| <qi+ql for £ < |u| < tt. 

Proof. It is noticed in [3J Lemma 2, Section 3, Page 182] that there exists a 
constant qi <1 such that |/(e*“)| < qi for e < |u| < tt. The lemma now follows 
from Corollary [S] and the definition of A (see (HOll b 

Due to m and Corollary [BJ we have the following lemma. 

Lemma 10. There exists positive constants s, 62 , 01,02 such that for |u| < e, 

|/t(e“)| < 1 - <52 v1 ^ + 


Lemma 11. For any r > 0, t —>■ + 00 , u —>■ 0 such that ut —>■ 0, 


At+r{en - Men = Men - /*+.(e“) = 


2rn{u)e 

(1 — g-ta(u)'j2 


(1 + ei(t,u)). 


where |ei(t,u)| < 0 (l)(ttt + + -v/ulnt). 

Proof. Due to Lemma (T] we have that ft{z) — ft+r{z) = At{z) — At+r{z). 
In order to prove the lemma, we use the following ingredients from the proof 
of Lemma [3 From equations ( 8 ),(9),( 10 ),( 12 ) of Lemma (3 we have At{z) = 
af{z)/if{t, z) + £ 5 ( 1 , z), where 




1 _£-*“(“) (l + 0 (l)ut) 
2 a{u) 


(1 + 0 ( 1 )i/m). 


t-i 

£ 5 ( 1 , z) = l/fiz) - ^ af'{z)Ak{z)Ok{l){z). 


(1 + O(l)Mt). 






From here we have that 


\s 5 {t + r,z) - e 5 {t,z)\< ^ \a^(z)Ak{z)0k{l){z)\<0{l) ^ ^ 

k—t 

\s 5 {t,z) < 0(1) Int. 

Since also \tp{t, z)\ > 0.5t(l — 0{t^/\u\)) > 2\e5{t, z)\, we have that 


t+r 


(z) 


t+r 


{z) 


< 


i>{t + r,z) + e 5 (t + r, z) i>{t + r,z)+ e^it, z) 

e“(*+’')“(“)(l + 0 (l)ut)|e 5 (t + r,z) - es,{t, z)\ 


< 


(|^/>(t+ r,z)| - \e^{t + r,z)\){\tl){t + r,z)\- \e^{t,z)\) 

\a^{u)e 


< 


<0{l) 




t|l-e-‘“(“)(l + 0(l)Mt) 


|2 ■ 


(42) 


Now let us consider the function 

a^{z) 2a{u)e~^°‘^'^\l + 0{l)ux) 

^ tli{x,z) + e^{t,z) 1 — e“““(“l(l + 0(l)ua;) + 0(l)a(u) Int 

One can see that x appears equally as a factor before real and imaginary part 
of 7 (a:). Since 7 (x) is also smooth for x > 0, for some x G it,t + r) we get that 


7 (t + r) — 7 (t) = r^'(x) = 2ra(u)e '““^“^(O(l)u — a(u))(l + O(l)ux)- 


[f _ g xa(u)(^i _j_ ( 9 (i)yi) -f- 0(l)a(u) Inf”*" 

g-xa(u) -2ra2(u)e-"^“(“) 

H-= - — - 

(l — e“““l“l(l + O(l)ut) + 0(l)a(it) Inf)^ (l — 


(l+ei(x,u)), 

(43) 


where |ei(x, u)| < 0(l)(ux+ + \/H Inx). The lemma now follows from (|43l) 

and (l42ll and the condition that t -+ +oo. 


We can get the following corollary from Lemma [TT] and Lemma |51 

Corollary 7. For any r > 0, t, N ^ +oo,u -+ 0 such that ut ^ 0 and N < 
2t\nn, we have 


ft^ienift+rien-Menf 


4x^0^(u)e 

(l_e-M.))4 


^—Na{u) 


1 + \ 
1 _ g-ta{u) j 


where |e 3 (t,u)| < 

















Proof. Due to Lemma [TT] and Lemma [51 it remains to bound the residual term 
as follows. 


^ut + + v^lnt^ + A^lnt ^|m| + < 

< 0(1) In^ t + ^/u + 2|u|t+—^ < 


< 


0(1) In^ t < 0(1) ^Inn ^\/|u|t + ■ 
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